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ABSTRACT 

We discuss the spin evolution of pulsars in the case where a superfluid component 



of the star is coupled to the observable crust on long, spindown timescales. The 
momentum transfer from the superfluid interior results in an apparent decay of 
the external torque and, after a dramatic increase, to an asymptotic decrease 
of the generic value of the braking index, e.g. n = 3, to values n ~ 2.5 if the 
magnetic field of the star does not decay over its lifetime. In the case where an 
exponential decay of the magnetic field towards a residual value occurs, the star 
undergoes a spin-up phase after which it could emerge in the millisecond sector 
qq ' of the P-P diagram. 
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The spin evolution of radio pulsars is governed by the braking torques acting on the 
solid component of the neutron star. A number of processes contribute to the energy losses, 
e.g. magnetic-dipole radiation (MDR) and a relativistic magnetohydrodynamic wind, which 
predict a braking torque oc —7 Q n , where n is the braking index, 7 - a positive constant, 
whose explicit form depends on the emission process assumed. Timing observations provide 
the values of n for four young pulsars (Kaspi et al. 1994; Boyd et al. 1995; Lyne, Pritchard, 
& Graham-Smith 1988; Lyne et al. 1996), which fall in the range 1.4 < n < 2.84, i.e. 
less than the canonical MDR value, n = 3. It is generally agreed that the observed 
grouping of pulsars in the P-P plane requires torque evolution in terms of variation of a 
pulsar's magnetic moment. The precise form of the time-dependence is, however, unsettled. 
The possibilities considered range from magnetic field decay on timescales of order 10 6 
yr (Lyne, Manchester, & Taylor 1985), decay on the same timescale towards a constant 
residual value several orders of magnitude lower than the initial one (Kulkarni 1986), to 
a marginal field decay (Romani 1990a). The separate class of millisecond pulsars, which 
have anomalously small braking torques, suggests that substantial torque decay occurs over 
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the evolutionary timescale. Their locations in the small P and small P part of the P-P 
diagram are commonly attributed to evolution in binary stellar systems, whether or not 
such pulsars are found in binaries today. A critical part of the evolution is the transfer of 
angular momentum and the attenuation of the magnetic field by accretion (Romani 1990b). 
Isolated millisecond pulsars may be accounted for by the evaporation of the companion by 
the pulsar's relativistic wind (Phinney et al. 1988; Ruderman, Shaham, & Tavani 1991) 
or three-body break-up processes (Romani, Kulkarni, & Blandford 1987). Detailed models 
suggest that evaporation may not operate quickly enough to account for all such objects. 
Though there is compelling evidence for the binary evolution scenario of millisecond pulsars, 
the existence of solitary objects among this class of pulsars appears somewhat enigmatic. 

The present paper studies, phenomenologically, the secular evolution of a solitary 
pulsar with a small fraction of interior superfluid coupled to the normal component on 
evolutionary timescales t s ~ P/2P ~ 10 3 — 10 6 yr. This assumption will be taken as a 
working hypothesis throughout most of the paper; the microscopic physics that might allow 
such a circumstance is discussed at the end of the paper. The implications of our hypothesis 
of weakly coupled superfluid interiors on the spin evolution in an accreting system are not 
discussed here. 

We shall adopt a non-linear version of the familiar two-component model of superfluid 
rotator (Baym et al. 1969). The equations are first derived for an arbitrary configuration 
in equilibrium and are specified for the spherical case thereafter. The Euler equations 
for the superfluid and normal components in the two-fluid superfluid hydrodynamic (e.g. 
Khalatnikov 1965) read 

Ps Vv? + f (m/) , (1) 

p„Vv?-f (m/) +& x \ (2) 

where v and p denote velocity and density; the subscripts s and n refer to the superfluid 
and the normal components; p is the pressure; </?, the gravitational potential, and ^ ex \ the 
external force acting on the normal component. It is assumed that the densities of the 
superfluid and normal components are conserved separately. Here the mutual friction force 
is written in general form as 

f(-/) = -{3[ ux [wx ( Vn -v s )]]-/?'[u,x (v n -v s )], (3) 

where u> = V x v s is the superfluid circulation; v = uj/uj] and the vortex tension forces 
(oc V x v) are dropped (we also ignore the effects related to the lattice stress energy; 
see e.g. Chandler & Baym 1986 for a discussion of these effects); the parameters (3,(3' 
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are phenomenological mutual friction coefficients. The macroscopic equations of motion, 
derived from (|1|) and (0), assuming hydrostatic equilibrium, are 

j t ( /( % - 4f) 0.i = K l + K>, (4) 
| (/(%• - iff) n nj = -Ki - K- + K\ ex \ (5) 

where Q s and Q n are the angular velocities of the superffuid and the normal components 
and iff 1 ** = J P( s /n)XiXj dV, are the moment of inertia tensors of these components; 
(/(»/») = Tr lff n) ). The internal torques K, K ', are defined as (summation over repeated 
indices is assumed) 

Ki = -2/3|O s | [^^- - 7^50^) + (4 S) ^^ fe - 4^(50^], (6) 

K\ = 2(3'\n s \u k e ijl I^8n j , (7) 

where SQ = fl n — £l s is the difference between the normal and superffuid angular velocity 
vectors. We suppose, further, that the neutron star decelerates under the action of magnetic 
dipole radiation, in which case n = 3 and 7 = B\ R 6 /6c 3 , where, B± is the surface magnetic 
induction component perpendicular to the pulsar's spin axis and R is the radius of the 
corotating star-magnetosphere system. Finite size corrections to the dipole radiation reduce 
the braking index to values between two and three (Melatos 1997), but as a fiducial value, 
we adopt the MDR value, n — 3. From now on we shall assume a spherical symmetry, 
in which case torque K vanishes. The resulting coupled equations can be integrated 
numerically for known initial conditions. Note that the linear model of Baym et al (1969) is 
recovered form these equations by defining a time-independent relaxation time (2|f2 s |/5) _1 . 

Shortly after birth, a cooling neutron star is expected to develop superffuid regions 
which carry most of the angular momentum of the interior in the form of an array of 
vortices, each carrying a single quantum of vorticity. When vortices nucleate, the superffuid 
and normal components are corotating, which is the initial condition for the subsequent 
integration of equations of motion for secular dynamics. As the external torque acts on the 
outer crust of the neutron star, the resultant spin down is communicated to different shells 
in the superffuid on different timescales determined by values of (3. Example evolutionary 
tracks of pulsars on the P-P diagram are shown in Fig. 1. 

Initially, the normal component spins down under the sole action of the external 
torque; the superffuid component remains uncoupled and, therefore, preserves its angular 
momentum acquired from the corotation phase at vortex nucleation. The imbalance 
between the superffuid and the normal component eventually starts to relax after the star 
evolves to times comparable to the dynamical coupling time ~ (f2o/5) _1 , where Qq is the 
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Fig. 1. — Evolutionary tracks of pulsars on the P-P diagram for magnetic field values of 
B ± = 10 13 and B ± = 2x 10 12 , a = 21^/31^ = 0.1 and log/3 = -14 (thinnest line), -16, -18 
(thickest line). The initial spin period is 10 msec. The labeling of the thick dots corresponds 
to log t (yr) with equal-time dots connected. 



- 5 - 



initial corotation frequency. This is marked by the "drop" of the evolutionary tracks in 
the P-P diagram: the decrease of the angular momentum of the crust reduces due to 
momentum imparted from the interior. For values log/5 > — 12 (small coupling timescales) 
the star effectively maintains corotation between all superfluid and normal components. 
The slope of the tracks is then determined by the external torque and is, therefore, exactly 
the same as it would be for the rigid body system. However, if there is weakly coupled 
superfluid with large coupling timescales, there is a memory effect in that the star reaches 
a given point in the diagram later than in the rigid body case. For very small (3's, the 
corotation phase is never achieved during the active lifetime of the pulsar, and the star 
crosses the deathline with P, P parameters radically different from the rigid-body caseQ. 

The time evolution of the apparent braking index is shown in Fig. 2. The deviations 
from the canonical rigid body, MDR value (n = 3) are nonmonotonic: the momentum 
impart from the superfluid causes a substantial increase of braking index because of the 
decrease in P(t) as a result of near cancellation of the external and internal torques on the 
dynamical coupling timescale. Subsequent evolution under the action of both internal and 
external torques implies braking indices in the range ~ 2.5 — 2.8, i.e. less than the canonical 
value 3 and in agreement with those deduced from the timing data; (these values are still 
considerably larger than the value n = 1.4 reported for the Vela pulsar (Lyne et al 1996); 
the case of the Vela pulsar is, presumably, exceptional because of the high glitch activity of 
this pulsar). The evolution of n(t) approaches the n = 3 value asymptotically. 

Relaxing the assumption of constant magnetic field, we next assume an exponential 
decay of the initial field to some residual value, i.e. we replace B± — > B± exp (—t/rd) + B™ s , 
where r d is the decay constant, B T f s (<£. B±) is the residual magnetic field value. The 
constant B± case is recovered in the limit — ► oo; in this respect a smooth crossover 
between two version of magnetic field evolution is expected as is varied. In particular, 
one would expect that the power low field decay case, implying a slower field decay rate, 
would lead to a result intermediate between those considered here with no new qualitative 
modifications. 

Simulations of the evolutionary tracks for = 10 6 yr and B± = 10 13 are shown in 
Fig. 3. For fixed and log/9 > —12 the evolution follows the rigid body case with an 
exponential decay on the timescale r^. For log/9 = —16 and —18 the slow coupling to the 
superfluid causes an initial drop of the tracks, which is followed by a zero crossing and a 



1 The theoretical deathline, beyond which the star turns off as radio pulsar, depends on the structure of 
surface magnetic and the radio-emission mechanism assumed, see K. Chen and M. Ruderman (1993); its 
precise location is, however, not needed for the present discussion. 
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Fig. 2. — The time evolution of the braking index for B± = 10 13 , a = 0.1; Solid, long- 
dashed, dashed and dotted curves correspond to the values log (3 = —12, —14, —16, —18. 
The initial spin period is 10 msec. 
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Fig. 3. — Evolutionary tracks of pulsars on the P-|-P| diagram with magnetic field decay for 
initial magnetic field B± = 10 13 and residual field B r f s = 10 10 G. An initial period 1 ms is 
assumed, while a is fixed at the value 0.1. The dashed and dashed-dotted curves correspond 
to rigid body evolution with and without magnetic field decay. The long dashed and solid 
curves correspond to tracks with log/3 = — 16 and —18 respectively. The arrows show the 
direction of motion (see text). The dotted curve is the fit to PSR 1937+21 with parameters 
given in the text. Equal-time thick dots labeled by logt (yr) are connected in some cases; 
whenever the direction of motion is ambiguous these are supplied by arrows. 
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spin-up; (the section labeled a on the example track with log/9 = —18 in Fig. 3). During the 
spin-up epoch (sections a upwards — b — c downwards) P < since the angular momentum 
transferred to the crust is greater than that lost via dipole radiation. The spin-up epoch 
is restricted by the residual field, which causes the tracks to cross zero once more (section 
c upwards), after which evolution enters a steady deceleration state (section d). When 
log/9 < —16 the initial drop immediately goes over to a spin-up stage, because the external 
torque has decayed substantially before the coupling between the superfluid and the normal 
component sets in, [i.e. (flfi)" 1 > r^)]. Note that the transitions from spin-down to spin-up 
and back to spin-down might require crossing the deathline: the pulsar will undergo 
period(s) of radio-silence, before eventually emerging as a millisecond pulsar. The possible 
set of parameters which fits the fastest known millisecond pulsar 1937+21 is not uniquely 
defined; the generic requirement, which does not seem to conform the present view on the 
interior of neutron stars, is the large moment of inertia of the weakly coupled superfluid 
(~ 60% of the total moment of inertia) and a fast field decay. In Fig. 3 the fit is obtained 
by setting the initial period to 0.8 ms, a — 1, log/5 = — 18, B± s = 10 9 and = 10 4 yr. 

To what extent can values of the phenomenological parameters a and f3 that produce 
spin up behavior be consistent with the microphysics of superfluid interiors? The a's 
are fixed by prescribing the equation of state and the composition of the matter along 
with the density-profiles of pairing gaps. The /3's, in the present model, are related to 
the local (microscopic) vortex-quasiparticle viscous friction coefficient rj, via the relation 
j3 = rip s u)/\rf + (p s uj) 2 ]. Of particular interest is the limit i]/p s uj — > 00, which shows that 
the strong coupling of quasiparticles to the vortex lattice on the local scale results in small 
values of f3, long coupling timescales, and effectively weak coupling of the superfluid to 
the normal component on macro-scales. The term "weakly coupled superfluid" is used 
throughout the paper in the latter sense. The opposite limit, p s ou/r) — > 00, which also leads 
to weak coupling, seems to be less plausible. 

Large friction and long coupling timescales in the superfluid crusts are predicted by 
the models of vortex creep (Link, Epstein & Baym 1993, Alpar et al 1996), where the 
coupling depends exponentially on the temperature. In this case the secular and thermal 
evolutions are coupled (e.g. Shibazaki and Mochizuki 1995), (3 is generally a function of 
time, and, therefore, our arguments apply only qualitatively. The minimal values of /3's, 
inferred from the respective relaxation times, lie in the range —12 > log f3 min > —14 for 
the time span 10 4 < t < 10 6 yr. For these values of (3 and typical values a < 0.03, the 
crust has a minor effect on the evolutionary tracks; a sizable effect, however, is expected for 
altering the apparent braking index (the value of a used in Fig. 2 is roughly by a factor 
3.3-10 larger than that expected for the crusts, but this can be compensated by choosing 
a slightly smaller initial spin period). A possibility, that would introduce a new timescale 
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in the problem, but is not incorporated in the present model, is the large scale motion of 
crustal matter and its superfTuid phases (Ruderman 1991). 

The dynamics of rotating superfluids in the quantum liquid core of a neutron star 
is well described in terms of the two-fluid hydrodynamics and eqs. (§]), (0) are strictly 
applicable. Among several frictional processes (Sauls, Stein & Serene 1982; Alpar, Langer, 
& Sauls 1984, Sedrakian et al 1995) the largest values are found for the process of electron 
scattering off proton vortex clusters forming a substructure of the neutron vortex lattice 
(Sedrakian et al 1995). The minimal value of (3 for this process is (3 m \ n ~ 10~ 18 at density 
p min ~ 8.3 x 10 14 g cm -3 , and corresponds to coupling times ~ (2Q (3)~ 1 = 2.5 x 10 6 yr for 
Qo = 6283 s _1 . This value places a lower theoretical constraint on the /3-parameter for the 
processes we have examined^. 

The BCS pairing would be suppressed when the velocity difference between the 
superfluid and the normal component exceeds the critical velocity v CI ~ (A/ep)vp, where A 
is the gap, ep and Vp are the Fermi energy and velocity. According to the familiar Landau 
criterion the quasiparticle excitations are spontaneously generated when the kinetic energy 
associated with the Galilean transformation to the moving frame exceeds the energy barrier 
~ 2A. Similar criterion for Cooper pair disintegration is obtained when considering the four 
point nucleon-nucleon scattering vertex for a pair moving relative to the background nuclear 
medium (Aim et al. 1996). The maximal relative velocity for the case with no magnetic 
field decay is <5t> max = <5fi max r ~ 6.3 x 10 8 cm/s, where r ~ 10 6 cm and 5Q max is the maximal 
angular velocity departure (corresponding to the case where superfluid rotates at the initial 
angular velocity while the normal component is virtually at rest). In the case of magnetic 
field decay, the maximal departure is determined by the field decay constant and we find 
that <5t> max < 10 7 cm/s in all cases. An order of magnitude estimate of the critical velocity 
at the nuclear saturation density n = 0.16 fm -3 gives 2 xlO 8 cm/s for A = 1 MeV. When 
■^max > Vcr the condensate will lose momentum via the pair-braking processes until the 
angular velocity departure drops below the critical value. Some models might maintain the 
situation where <5t> max < v CT which would imply earlier and smaller "drops" of evolutionary 
tracks. 

Observationally, our conjecture may be testable. First, if the currently measured 
braking index is significantly influenced by internal torques, then it should change over time 
in a way that can be tested statistically in a sample of pulsars (Cordes & Chernoff 1997). 



2 For the neutron star model with canonical mass 1.4M based on the AV14 equation of state of Wiringa, 
Fiks and Fabrocini (1988) we find a — 0.094, for the case when superfluidity persists through the core density 
region p m ; n < p < p c , where p c — 1.03 x 10 15 g cm~ 3 is the central density. 
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Second, the evolution towards the sector of millisecond pulsars in the P-P diagram does 
not imply increase of the mass of pulsar, while the pulsars spun-up in an accretion phase 
should be more massive than their non-spun-up counterparts by ~ O.1M . The duration 
of the spin-up phase in the present model for internal torques is the superffuid relaxation 
time, which e.g. for times ~ 10 6 yr implies that in a population of Gyr old millisecond 
pulsars there is a probability of < 10~ 3 of observing a pulsar with P < 0. 

To summarize, we have shown that weakly coupled superfluid in pulsars affects the 
secular dynamics of pulsars in two ways: (i) For constant magnetic field the internal torque 
causes an apparent degradation of the braking torque, thus mimicing field decay. The 
braking index, after the epoch of drastic enhancement, reduces to values ~ 2.5, close to 
those measured for young pulsars, (ii) In the case where the magnetic field is attenuated 
the star undergoes a spin-up phase after which it emerges as a solitary millisecond pulsar; 
one should note, however, that fits to the fastest known solitary millisecond pulsars require 
unusual parameter values which seem to be incompatible with the present view on the 
structure and field evolution of neutron stars. 
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